Angle resolved Cu and O photoemission intensities in Cu02 planes 
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Using a mapping from the three-band extended Hubbard model for cuprate superconductors into 
a generalized t — J model, and exact diagonalization of the latter in a 4x4 cluster, we determine 
the quasiparticle weight for destruction of a Cu or an O electrons with definite wave vector k. We 
also derive an approximate but accurate analytical expression which relates the O intensity with 
the quasiparticle weight in the generalized t — J model. The fc-dependence of Cu and O intensities 
is markedly diflferent. In particular the O intensity vanishes for k — (0,0). Our results are relevant 
for the interpretation of angle-resolved photoemission experiments. 



I. INTRODUCTION 

While the wave vector dependence of the quasiparticle 
weight in quantum antiferromagnets has been studied in- 
tensively since the discovery of high- Tc superconductiv- 
ity, the interest on the subject has been revived by the 
angle-resolved photoenaission (ARPES) experiments on 
insulating Sr2Cu02Cl2tl. Several theoretical works ap- 
pear|ecLfitting the observed dispersion using generalised 
i — Jd^eI, spin-fermioncl or one-band Hubbard modelsEJO. 

More rec^r^bi, ^e photoemission intensities have been 
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discusseS^ 
thfirt — J mode 

and Sushkov et 



i, compared with previous results in 
^' ^XL, particular Lema and one of 
have developed two differ- 



ent methods to calculate the quasiparticle weight in the 
generalized t — J ot one-band Hubbard models in the 
strong coupling limit, using the self-consistent Born ap- 
proximation (SCBA). The results of both methQdjS, an 
analytical approach based on the "string. picture" E3, and 
exact diagonalization of a 32-site clusteriO were compared 
recentlyls. The method of Sushkov et al. introduces 
spurious low-energy peaks in the Green function which 
can, however, be identified and eliminated. The other 
SCBA method compares better with exact diagonaliza- 
tion and the results of the string picture underestimate 
the weights. However, since the operators c^^ entering 
generalized t — J ov one-band Hubbard models are effec- 
tive operators which cannot be trivially translated into 
Cu and O holes of the original system, the above men- 
tioned efforts are insufficient for a comparison with ex- 
periment. 

Experimental evidence 
in cuprate superconductorst 
density functional calculatj. 

band Hubbard model H^,]^'^^ as the starting point 
for the description of these systems. iJsf, contains_Cu. 
3da;2_j,2 and O 2p(jj:irbitals. To explain some RamarSa 
and photoemissiono experiments at excitation energies 
aboverieV, it is necessary to include other orbitals in the 
modelEa, but these are not important for the energy scale 
of the ARPES experiments of Ref. 1. However, even re- 
stricting the basis to the above mentioned orbitals, the 



the symmetry of holes 
s, well as constrained- 
JisiPq justify the three- 




size of the systems which can be diagonahzed numericajlly 
at present, do not contain more than four unit cellsEH. 
This is too sipall to discuss the above mentioned ARPES 
experimental. Thus, to study this problem by numeri- 
cal methods at zero temperature, it is necessary to inte- 
grate out the high-energy degrees of freedom. Further- 
more, analytical approximations like slave bosons give 
better results_svhen applied to an appropriate low-energy 
HamiltonianEfl, and the successful SCBA cannot be ap- 
plied to H^b. 

Several— Jossf-energy reduction procedures have been 
proposedEa~LJ. Eliminating the Cu-0 hopping tpd by 
means of a canonical transformation, leads tp-the spin- 
fermion Hsf (or Kondo-Heisenberg) modcSEJ. Al- 
though tpd is in principle not small enough to guarantee 
the accuracy of the resulting i?s/, this effective Hamil- 
tonian, with parameters renormalized to fit the energy 
levels of a CUO4 cluster, reproduces very well optical 
and magnetic properties of H^b in a—GujjOg clusterE3. 
Also 0|iiepb|and generalized HubbardoO'Ej and t — J 
modelsLj'O were derived. The latter represent the high- 
est low-energy reduction reached so far, and after the first 
proposal of Zhang and RiceES, further work confirmed 
that a generalized t — J model Hqij reproduces.. accu.- 
rately the low-energy physics of the other modelsE^rLjEJ. 
In particular, projecting the Hilbert space of Hsf onto 
local (non-orthogonal) Zhang-Rice statescll, mapping the 
model in this reduced Hilbert space to Hctj, and solv- 
ing numerically the latter, we obtained a band structure 
and magnetic properties which agree very well witiL.the 
corresponding properties calculated directly on HsfLj. It 
is important to emphasize that to calculate any property 
of the cuprates, expressed in terms of the expectation 
value of an operator of H^b, using an effective low-energy 
Hamiltonian, the mapping procedure should be extended 
to the operator of the quantity to be calculatedLHlHcJ, or 
alternatively the relevant states of the effective Hamil- 
tonian should be mapped back onto the corresponding 
states of H^b- 

In this work we calculate the low-energy part of the Cu 
and O ARPES, using the low-energy reduction from H^b 
to Hsf and from it to Hqij, mapping the local Zhang- 
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Rice singlets to vacant siteai^'cS. The relevant operators 
of H^h are mapped to -ffs/, and the ground state of H^f 
is constructed from that of Hctj in a system contain- 
ing 4x4 unit cells. For the O ARPES we give a simple 
recipe to relate it with the quasiparticle weight in Hatj, 
which can be calculated with the SCBA or other analyt- 
ical approachea3. We find significant differences between 
Cu and O ARPES. Since the photoeniission cross-section 
for Cu d and O p orbitals have different dependences 
on the incident energy of the photorO, these differences 
should be accessible to experiments. 

In section II we briefly review the mapping procedure 
and derive the equations necessary to express the ARPES 
results in terms of numerical or (in some cases) analytical 
results on Hctj- Section III contains the results and 
section IV the conclusions. 



II. MAPPING PROCEDURE AND RELEVANT 
EQUATIONS 

A. Deriving the spin-fermion model from H^b 

Our starting point for the description of the supercon- 
ducting cuprates below leV is the extended three-band 
Hubbard model. To simplify the writing we change by 
-1 the phases of half of the O and Cu orbitals in such 
a way that the hopping matrix elements do not depend 
on direction. The original phases should be restoredrin 
the comparison with the experimental ARPES resultsEJ. 
The Hamiltonian takes the form: 

icr ja 

i j 

+ Upd ^ d\ad^ap\+f,„,P^+5a' 
iSaa' 

+tpd [p]+sadi<y + h.c. 

16 a 

-ipp^Pt+taPj^ (1) 

where d\^ {p\cr) creates a hole on the Cu 3d (O 2p) orbital 
at site i (j). The four nearest-neighbor (next nearest- 
neighbor) O sites to Cu site i are denoted hy i + 5 (i-t-7). 
A canonical transformation which eliminates terms linear 
in tpd, retaining alsp-Jjiie fourth order terms, leads to the 
spin-fermion modelME^: 

Hsf = Yl pl+S'aP^+S'y {{ti + t2)d + 2S. • S.+i) - t2] 

~^PP X! + ^(Sj • - |) 

j^fo- id 



+ ^^(S,.S„+2.-- J), (2) 

iS 

where P^^ are effective O creation operators, and Si 
{^i+s) is the effective spin at Cu site i (O site i + 5). 
Due to the fact that tpd is not very small compared to 
A or J/^i — A, the expressions for the parameters of Hsf 
obtained from the canonical transformation up to fourth 
order in tpd are not accurate enough. However, this short- 
coming is avoided if the parameters of H^f are renormal- 
ized to fit the energy levels of H^i, which in the limit 
tpd corresponds to a level of Hsf, in a CUO4 clus- 
ter with one and two holes. Since the case tpp ^ has 
not been described before and the information is neces- 
sary for the expressions of the ARPES results, we briefly 
review this method. 

For two holes in the CUO4 cluster the 16 eigenstates of 
Hsf can be classified in four spin singlets and four spin 
triplets distributed in six energy levels: one Pi (invariant 
under the point group operations), one (transforming 
like — y^) and a doublet (transforming like x, y) for 
each total spin. The spin multiplicity 2S-I-1 will be de- 
noted in the superscript. The ground state is the invari- 
ant singlet (TI)) = ^E^sipl+sAl - pI+sAM' 
which represents a Zhang-Rice singlet. Each eigenstate 
of Hsf has a corresponding eigenstate of Hsb, which is 
the lowest eigenstate of a small matrix in the correspond- 
ing symmetry sector T^. The largest matrix corresponds 
to P} and is reproduced here for future use: 



/ —2tpp 2tpd V^tpd V^tpd V&tpd ^ 

2tpd A — 2\/2tpp —2^/2tpp 

V2tpd -2V2tpp AO 

V2tpd -2y/2tpp A + f/p 

V%/8tpd Ud-A-2UpdJ 

(3) 



The basis states of Eq.(|^) are the following: 

1 2) = ^Y.^pI+sM+RH- pI+sipI+b.5^M^ 

|3> = ^EpWUio>> 

|5)=4t4|0): (4) 

where R6 is the result of rotating 6 by tt/2. To obtain 
the optimum parameters of Hsf., we adjust them to fit 
exactly the three lowest energy levels (P}, Pg and Pf), 
the highest one (Pf), and the average of the other two 
((P^ + Tl)/2). Calling Eat the lowest energy of Hsb in 
each symmetry sector, the result is: 
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PP 4 

E3biTl)~E,,{rl) + 2t'p^ 

^ 

^ ^ 8 2 4 

Jk = 2iti + t2) + E3biTl) - Esbirl) (5) 

As an example for the parameters of H-^b for La2Cu04, 
obtained frcm constrained-density-functional approxima- 
tion in RefEl {Ud = 10.5, Up = 4.0, Upd = 1.2, A = 
3.6, tpd = 1.3 and tpp = 0.6, all energies in eV) , we 
obtain t'^^ = 0.56, h = 0.37, t2 = 0.08, Jk = 0.62. The 
value oO, which is affected by other orbitals noitincluded 
in iJsbCEl is taken as J = 0.13 from experimentcZl. 



where \a'\ < 1, because part of the spectral weight of 
pI^ is distributed in high-energy states which are out of 
the Hilbert space of Hsf. Eqs.(||) and (^) were shown 
to be accurate anough in previous comparison of the Cu 
and O rphotoemission spectra of H^b and Hgf in a CU4O8 
clusterE3. Here, to calculate a', we solve exactly H^b and 
Hsf in a CU2O cluster including and O atom and its two 
nearest- neighbor Cu atoms, with two and three holes. 
The first (second) member of Eq. (H) is applied to the 
S=0 ground state of H^b {Hsf) with two holes. For Hsf, 
the result is a linear combination of two eigenstates with 
total spin S=l/2, which correspond to the low-energy 
part of the result for H^b- Then, a' is determined fitting 
the coefficients of these two states. Inside the range of 
reasonable parameters of H^b, we obtain |a'|^ ~ 0.44. 



B. Mapping of the operators 

We have to express the hole creation operators dl^ 
and pI^ in the basis of Hsf in order to calculate pho- 
toemission properties. In the lowest non-trivial order in 
the canonical transformation which eliminates tpd, one 
obtains for the dj-j. operator transformed into the spin- 
fermion basiS: 



(6) 

and similarly interchanging spin up and spin down, with 

nia — dl^dia- The values of a, 6, c which are ob- 
tained from the canonical transformation are not accu- 
rate enough. To improve them, we ask that all matrix 
elements of the second member of Eq. (||) between states 
of Hsf in a CUO4 cluster with one and two holes, should 

coincide with the matrix elements of d\^ between the cor- 
responding states in H^b- The result is: 

a = vl2, h = -u\A^\/y/?> + {l~\Ai\)vlA 
c = u\A^\ / \/% + {I + \Ai\) V / A, with u, u > 



u = — 
2 



A + iJpd - 2t 



PP 



2J{A + Upd-2tpp)' 



ml 



(7) 



The Ai are the coefficients of the ground state of the ma- 
trix Eq.(||) in terms of the basis set Eq. (0): jgsf, (T})) = 

i 

In the lowest non-trivial order in the canonical trans- 
formation, the transformed operator of p\^ is not 
changed. Following a similar procedure as above, we as- 
sume that it is a good approximation to use: 



Pi, 



(8) 



C. From Hsf to a generalized t — J model 

There is numerical evidencJl^ that in the low-energy 
eigenstates of Hsf , the O holes are iiV-the ground state of 
a CUO4 cluster (a Zhang-Rice singletE^I) . Defined in this 
way, Zhang- Rice singlets centered in nearest-neighbor Cu 
orbitals are non orthogonalo Using a jppjector P2 over 
these non-orthogonal Zhang- Rice statescil, P2HsfP2 can 
be mapped into a generalized t — J model Haj, in which 
each Zhang-Rice singlet at a CUO4 cluster, is replaced by 
the vacuum (no holes) in the cluster. Retaining the most 
important terms, Hctj takes the formeil: 



+*3 5Z ^1+2 Act 



17 CT 



iAct 

+t" E rfl+A'CT'^^^+ACT (1 - 2S. • S,+a) 
J 

+ 2 



J2 (s« • S,+A 

iA<T ^ ^ 



(9) 



where A = 25 (F = 27) are vectors connecting first (sec- 
ond) nearest-neighbor Cu atoms, and: 

t[ = (lO-it'pp + 246ti + 410t2 + 51Ja') /512 
4-(l3t;p-llti)/64 



t' 



-lKi/128 



t" = {8tpp - 18ti - 6t2 + 3Jfc) /256 



(10) 



As in Hsf there is an implicit constrain of forbidden 
double occupancy at any site. For the typical param- 
eters of ifsf, mentioned above, Eqs. (^0|) give: t\ = 0.42, 
= 0.05, t'^ = 0.06, t" = 0.003. 

The low-energy eigenstates l^'^j) of Hsf can be ob- 
tained from those l^^t j) HctJ simply by dressing the 
vacant sites with Zhang-Rice singlets: 
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T\^htj) 



\T^T\ 



where Ni. — 



E (^!+^T 4- p^+si 4) (1 - + 



(11; 



This equation, together with Eqs.(g) to (||) allow to cal- 
culate the Cu and O photoemission spectra of the three- 
band Hubbard model H^^ from the eigenstates of the cor- 
responding generalized t — J model Hqij. For the sake 
of clarity, and in absence of a more detailed knowledge 
about the experimental situation, we neglect effects of 
interference and the dependence on the polarization of 
the incident radiation and direction of the photoemitted 
electron. Then, in the insulating state, the Cu and O 
contributions to the intensity of the lowest ARPES peak 
are given by (the contributions for both spins are the 
same) : 



and in the sun 
over j runs over half the O atoms (those which contains 
2pa; orbitals or in other words, their nearest-neighbor Cu 
atoms lie in the x direction). The norm ^ 1, due to 
the non-orthogonality of Zhang-Rice singlets centered in 
nearest- neighbor Cu sitescil. Since l^'^jj) contains only 
one vacant site then (1 — ni){l — nj)\^Q^j) 
Using this and Eq.(|ll|) one has: 



for i ^ j. 



'ctj) 



Io{k) = 2 



di 



kCT 



Pkxa 



*0 



= ^MUi + ^~UJI*'*^) (15) 

where the last equality makes use of the fact that rii = 

or 1 and E (1 - "0 = 1- From Eqs.(|l|) to (|l|) one 

■i 

obtains: 



.,2 f 



Pkya 



(12) 



Im (k) = 



2Nk 



-Zi (-k) 



(16) 



where |*o) is the ground state of H,f and Hgu in the since from symmetry (k)^ (~k) = Z^ (-k) , 



insulating system, l^'^y) the lowest energy eigenstate of 
Hgf for one added hole (which leads to a non-zero ma- 
trix element) and dj^.^, pl^^^, and pl_y^, are the Fourier 
transforms of the three creation operators of a unit cell. 



D. O intensity vs quasiparticle weight in Hctj 

The formalism presented in the rest of this section al- 
lows us to calculate, in the next section, the Cu and O 
contributions to ARPES, from the eigenstates obtained 
from exact diagonalization of finite systems. To calcu- 
late the Cu part with analytical approximations applied 
to Hqij requires further algebraic elaboration which is 
beyond the scope of this work. However, as we show 
below, there is a simple analytical relation between the 
O contribution, generally the most important, and the 
quasiparticle weight of Hctj for one added hole. The 
latter jquantity has been calculated accurately with the 
SCBAQ'Erta and comparejd with results of other analytical 
and numerical methodsEl. 

The quasiparticle weight in Hctj is: 



we have for the total O intcnsityc3: 



/n(k) 



cos^ (^) 



C0S2 ( ^ 



V|2 



-^T (k) (17) 



^.(k) 



*0> 



(13) 



while the contribution to the intensity from, for example, 
2pa; orbitals and spin up is (Eqs.(||),(|ll]), and (p^): 



JS,(k) =|af |(*o|Pk.Tl* 



In the cluster of 4x4 unit cells, we obtain iVk = 0.36 for 
all k. with error less than 10%. The small dependence 
of the norm on wave vector is to be expected 

in an antiferromagnetic background for realistic parame- 
ters of HctJ- rAs it becomes particularly clear within the 
string pictural3, the motion of a vacant site in a quan- 
tum antiferromagnet can be divided into a fast motion 
around a fixed position on the lattice, on the scale of ~ 3t, 
against a string linear potential created by the distortion 
of the antiferromagnetic order, and a slow motion of the 
polaronic cloud, which determines the quasiparticle dis- 
persion (with a width ^ 2 J). N-^ is clearly determined 
by the physics inside the polaronic cloud and is thus es- 
sentially independent of its wave vector k. 

From the above discussion, it is clear that the wave 
vector dependence of Iq (k) is given by that of the 
quasiparticle weight of Hctj, and the factor cos^ (■^) -|- 

cos^ ^-^^0. This factor is very important and leads to 

the fact that for wave vector (tt, tt) in the notation of 
Eq. (|l|) (k = (0, 0) when thCpOriginal phases are restored 
to compare with experimentllj), there is no O contribu- 
tion to the low-energy ARPES. This is particularly clear 
when the on-site O repulsion Up = 0. In this case, from 



1 ^' 



Vlv 



Eq. (§, 



t 

(7r,7r)a(T 



(7r,7r)c 



with a — 



X or y, 



Nk 



(14) 



I.e. vl s 



does not hybridize with the Cu 3a;2_j,2 or- 



bitals. Then all the O weight resides in a well defined 
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quasiparticle at energy A 3.6 eV, while the low-energy 
quasiparticles, involved in the formation of Zhang-Rice 
singlets, lie at negative energies (with the zero of one- 
particle energies of Eq. (tthEj) . 



III. RESULTS 

In this section we present the result of exact diago- 
nalization of Hqij as an effective model representing the 
low-energy physics of H^f,, in a system containing 4x4 
unit cells. At the end we use Eq.(|l^ and previous re- 
sults of Z| (k) to obtain Iq (k) in larger clusters. 
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FIG. 1. Quasiparticle energies (top), oxygen intensity 
(middle) and Cu intensity (bottom) as a function efi wave 
vector, for parameters of Ho^i, calculated for La2Cu04Ej. The 
square symbols and error bars a|fc the top correspond to the 
observed ARPES in Sr2Cu02Cl2B. 

At the top of Fig. 1 we show the quasipacticle dis- 
persion Ak with the original phases restoredcj and in 
the electron representation (upside down with respect 
of the hole representation of H^b), to facilitate compari- 
son with experiment. The parameters of H^b were taken 
from Ref.EEl and those of Hctj were determined from 
the mapping procedure, except for J = 0.13 which was 
taken from comparison with Raman experimentscZI. Tak- 
ing into account that there are no fitting parameters, the 
agreement with the experimentaly measured dispersion 
in Sr2Cu02Cl2 is very good. The discrepancies around 
(tt, 0) can he, ascribed to some finite-size effects in the 
4x4 clustedlJ, and to the fact that the parameters of 
ffsb for La2Cu04E£l should differ somewhat from the cor- 
responding ones for Sr2Cu02Cl2. A consequence of the 
upward shift in Ak for k = (tt, |) ,(7r,0) and (f,0) is 



that the quasiparticle wei^t (k) of Hqij is exagger- 
ated for these wave vectorsB. This is due to the fact that 
for larger binding energy of the added hole, less magnons 
are excited and the quasiparticle is more similar to the 
bare hole, increasing Z(j (k). 

The O and Cu intensities given by Eqs.([l^) are com- 
pared in Fig. 1. As explained above, both intensities are 
exaggerated for wave vectors (tt, f ) jfvr, 0) and (f,0) . 
For lo (k) this is clear when Eq.(^@ with \a'\^ l^^vm 
1.22, and the weights of Hctj calculated by the SCBAO-eI 
are used. However, these finite-size effects do not af- 
fect the characteristic strong variation of the O intensity 
around the E line (joining (0,0) with (7r,7r)). Iq (k) is 
maximum for k = ? f ) f^nd very-small for k ~ (tt, tt) , 
as for the generalized t~ J modelETEHlJ. However, in con- 
trast to Za (k) for Hctj, Iq (k) vanishes at k = (0,0) . 
This is a consequence of the different symmetry of the Pa 
and da;2_y2 orbitals (or Zhang- Rice excitations) at that 
point, as explained (in different terms) at the end of the 
previous section. 

In contrast to Iq (k) , the Cu intensities for k = (f , f ) 
and k = (0, 0) are similar and rather large in comparison 
with other wave vectors. Since Iq (0,0) = 0, the exper- 
imental ARPES intensity at k = (0, 0) is determined by 
the Cu part. The maximun of /cm (k) for k = (f;0) 
is probably not realistic for the parameters of H^^ which 
correspond to Sr2Cu02Cl2, and should be reduced as the 
corresponding A^ o approaches the observed quasiparti- 
cle energy. 




K=(7t,7t) 




K=(7t,0) 



K=(7t,7i:/2) 



K= (71,71) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 

too (eV) 



FIG. 2. Oxygen (top) and Cu (bottom) intensities for sev- 
eral wave vectors as a function of 0-0 hopping. 

In Fig. 2 we show the evolution of several intensities as 
a function of the 0-0 hopping tpp of H^b- The importance 
of this term is that as tpp increases, the three site term 
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t" becomes positive (see Eg. (fi^ ) , particularly if i" is ob- 
tained by fitting energy levelsEB instead of the analytical 
expression Eq. (|l^) we used here. In turn, moderate 
positive values of t" favor a resonance-valence-bond su- 
perconductiug. ground state with (predominantly) da;2_j,2 
symmetryEJ'EJ. The effect of tpp on the intensities is to re- 
duce lo (k) and Icu (k) for k = (tt, f ) , (tt, 0) and (f , O) . 
Also Icu (0, 0) decreases with tpp. This is mainly a con- 
sequence of a shift downwards of the corresponding Ak- 
As a consequence, the dispersion, and also apparently 
the intensities, compare better with the ARPES results 
in Sr2Cu02Cl2, if tpp ^ 0.6 eV or larger. 
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FIG. 3. Quasiparticle dispersion (top) and oxygen intensity 
(bottom) as a functirai of wave vector using exact |diagonal- 
ization (open circles)!] and SCBA (solid triangles)El results. 
Parameters are J — —t'2 — O.St'i, tg = 0.2t'i and t" = 0. 

For the m incident energy used in the ARPES 
experimental! in Sr2Cu02Cl2 , the cross section for pho- 
toemiting Oh^p electrons is near two times that of Cu 
3d electronsE3. This fact and our previous results sug- 
gest that the observed intensity is given essentially by 
lo (k), except for k near (0,0). We have used Eq.(p^) 
with Nk — 0.36 constant, in order to relate Iq (k) with 
previous accurate results for (k) in Hctji-, exact di- 
agonalization of a square jpluster of 32 sitesQ, and the 
SCBA in a 16>t|16 cluster^!. The parameters of Hctj, 
taken from Ref.L!, are near the optimum ones for fitting 
the dispersion relation Ak, with t" ~ 0. The resulting 
Ak and Iq (k) are shown iu-jFig. 3. Due to the fac- 
tor sin^ {kx/'2) + sin^ {ky/2ytB, the intensities along the 
E line are asymmetric and smaller near the Brillouin 
zone center, contrary to what was observed experimcn- 
taly. We should state that a small admixture of the Cu 
Sd^*^ configuration in the ground state of the undoped 
system (of order t'i^/[{A + Upd)Ud], which we have disre- 



garded here) has the effect of increasing the Cu ARPES 
near k = (0, 0), as it is clear in the stro«g..jOenpling limit 
of the one-band Hubbard model (i?ib)l3ll2!'E^. However, 
clearly this effect is negligible on the O ARPES. The 
above mentioned asymmetry is even enlarged if a neg- 
ative t" = —J/4 (as that which com^form a canoni- 
cal transformation of Hit,) is includedQij. This suggests 
again that values of tpp ^ 0.6eV or larger, leading to 
positive t" , are more realistic. 

IV. CONCLUSIONS 

We have developed a formalism which allows to calcu- 
late separately the low-energy part of the angle-resolved 
photoemission intensity from either O 2p^ or Cu 3(1^2 _y2 
orbitals, using a generalized < — J model as an effective 
model for cuprate superconductors. This low-energy re- 
duction is the only way to calculate the wave vector de- 
pendence of the intensities by exact diagonalization of 
finite systems, since at present it is not possible to diago- 
nalize directly the three-band Hamiltonian in a periodic 
system large enough to contain the minimum necessary 
sampling of the Brillouin zone. 

For the insulating system, the O intensity can be very 
well approximated as 

lo (k) ^ 1.22 (k) (sin^ {kj2) + sm\ky/2)) 

where Zcr (k) is the quasiparticle weight of the effective 
generalized t~ J model. Thus, Iq vanishes at the F point 
k = (0, 0) . Since this is a consequence of the different 
symmetry of O 2pa- states and low-energy excitations at 
that point, this result should persist with doping. 

Our numerical results in a cluster of 4x4 unit cells, for 
parameters calculated for La2Cu04, show that Iq (k) is 
largest for k = ^) , and at that point, the Cu inten- 
sity Icu (k) is nearly three times smaller. Instead while 
/cu (k) has similar values at k — and near the 

F point, lo (0,0) = 0. The fact that lo (k) and Icu (k) 
dominate in different regions of the Brillouin zone, makes 
it possible to separate both contributions experimentally. 
For an analysis-jof the experiments, as those carried out 
in Sr2Cu02Cl2EJ, the separation in Cu and O contribu- 
tions is important, since the cross section for photoemit- 
ing electrons in O 2p or Cu Sd^^-y^ orbitals are differeajb 
and have different dependence on the incident energyM. 
For a quantitative comparison with experiment, it is nec- 
essary add the amplitudes (instead of the intensities) of 
the scattered waves from the three atoms per unit cell, 
multiplied by their respective scattering amplitudes, tak- 
ing into account the polarization of the incident photons, 
and the direction of the photoemited electrons. This does 
not require an extension of our formalism. In addition, 
for any particular scattering amplitudes and polarization, 
the expected trends can be extracted from the present re- 
sults. 
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Agreement with the observed intensities seems to im- 
prove for tpp > 0.6eV, which in turn faywijs-|an RVB 
ground state and d-wave superconductivitjC3'Ell. 
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